LOG-MAJORIZATION OF THE MODULI OF THE EIGENVALUES 
OF A MATRIX POLYNOMIAL BY TROPICAL ROOTS 
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Abstract. Wc show that the sequence of inoduh of the eigenvalues of a matrix polynomial is log- 
majorized, up to universal constants, by a sequence of "tropical roots" depending only on the norms 
of the matrix coefficients. These tropical roots are the non-differentiability points of an auxiliary 
tropical polynomial, or equivalently, the opposites of the slopes of its Newton polygon. This extends 
to the case of matrix polynomials some bounds obtained by Hadamard, Ostrowski and Polya for the 

C^ roots of scalar polynomials. We also obtain new bounds in the scalar case, which are accurate for 

^~H "fewnomials" or when the tropical roots are well separated. 
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1. Introduction. Let p{x) — X]?=o'^J^"'' "^i G C be a polynomial of degree 
n in a complex variable x. Let Ci, . . . , C„ denote the roots of p{x) arranged by non- 
I— I decreasing modulus (i.e., |Ci| ^ . . . ^ ICnD- We associate to p the tropical polynomial 

Ph tp{x), defined for all nonnegative numbers x by 






The tropical roots of tp, Ui, . . . , «„, ordered by non-decreasing value (i.e., ai ^ . . . ^ 

I " I Q:„), are defined as the non-differentiability points of the function tp, counted with 

certain multiplicities. They coincide with the exponential of the opposite of the 
'""' slopes of the edges of a Newton polygon, defined by Hadamard |Had93| and Os- 
jj^ trowski |Ost40a[ IOst40b| as the upper boundary of the convex hull of the set of points 
vQ {{j, log \aj\) I ^ J ^ n}. The logarithms of these roots were called the inclinaisons 
^s^ numcriques by Ostrowski. One interest of these roots is that they can be easily corn- 
ed puted (linear number of arithmetic operations and comparisons). See Section^ below 
.^ for details. 

f^ Hadamard was probably the first to prove a log-majorization type inequality for 

C^ the modulus of the roots of a scalar polynomial by using what we call today the 

^~~' tropical roots. His result (page 201 of |Had93| . third inequality) can be restated as 

^ follows in tropical terms: 



ai • • • ckfe fc + 1 

This bound, proved in passing in a memoir devoted to the Riemann zeta function, 
remained apparently not so well known. In particular, the special case |Ci|/ai ^ 1/2 
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is equivalent to the homogeneous form of the classical bound of Cauchy, established 
later on by Fujiwara |Fujl6| , and a weaker inequality, with a^ at the denominator 
instead of ai • • • a/j, appeared later on in the work of Specht [SpeSS]. 

Ostrowski proved several bounds on the roots of a polynomial in his work on 
the method of Graeffe |Ost40al IOst40b| , in which he used again the Newton poly- 
gon considered by Hadamard. In particular, he obtained the following upper bound 
(see [Ost40al §7]), 



which can be thought of as a generalization of a "reverse" of the Cauchy inequality 



due to BirkhofF |Birl5j (corresponding to the case fc = 1 in (1.2)). He also gave a 



different proof of a variant of (1.1), with the constant l/(2fc) instead of l/{k + 1), and 



reported a private communication of Polya, leading to a tighter constant 

IClC2- --Cfcl ^ 1 Q gN 

ai---ak ^ ^£(fc)(fc + l) ' 
with 

f(fc):-(^)'<e. (1.4) 

In this paper, we generalize the bounds of Hadamard, Ostrowski, and Polya, to 
the case of a matrix polynomial 

P{X)^Ao + AiX + --- + Ad\'^, A, eC"^" j = 0,...d. 

We now associate to the matrix polynomial P the tropical polynomial 

tp(x) := max IMJIx"' , (1.5) 

where || • || is a norm on the space of matrices, and show that the moduli of the roots 
Ci: • ■ • I Cnd of P can Still be controlled in terms of the tropical roots ai, . . . , q;„ of tp. 
Our results give bounds on the ratio |Ci ■ • ■ C,nk\/{oLi . . . a/c)", which extend and refine 



the above bounds. In particular, in Theorem 2.4 we extend the lower bound (1.31 



of Polya to the matrix polynomial case, and in Theorem |2.12[ we extend the upper 



bound ( 1.2 ) of Ostrowski. Moreover, we obtain other lower bounds that are new even 
in the case of scalar polynomials. In particular, in Theorem |2.H we obtain a lower 
bound which may be tighter for "fewnomials" . In Theorems |2.5| and |2.6[ we obtain 
general lower bounds, which extend the bound of Polya and its extension to the matrix 
case, and which may be much tighter when the tropical roots are sufficiently separated. 
Then, all together our results show that the tropical roots give tight estimates of the 
moduli of the eigenvalues if the tropical roots are sufficiently separated and if certain 
matrices are sufficiently well conditioned. 

The results of the present paper combine ideas from max-plus algebra and trop- 
ical geometry, and numerical linear algebra. In |ABG98[ IABG04[ lABGOSj , Akian, 
Bapat, and Gaubert studied the eigenvalues and eigenvectors of matrices and ma- 
trix polynomials which entries (the Aj) are functions, for instance Puiseux series, of a 
(perturbation) parameter. It is shown there that the leading exponents of the Puiseux 
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series representing the different eigenvalues (resp. eigenvectors) coincide, under some 
genericity conditions, witli certain "tropical eigenvalues" (resp. eigenvectors). This 
can be interpreted in the light of tropical geometry (see [VirOll IIMS07[ IRGSTOS] for 
introductions) , where one defines the amoeba of an algebraic variety V C K" over a 
field K, equipped with a valuation v, as the image of V by the map which applies the 
valuation v entrywise. When v is non-archimedean, in particular for the usual valu- 
ation (leading exponent) on the field of Puiseux series, Kapranov theorem [EKLOB] 
characterizes the amoeba of a hypersurface g = as the set of "roots" (nondiffer- 
entiability points) of the tropical polynomial function with coefficients equal to the 
valuations of the coefficients of q. Some of the results of ABG98, ABG04, A BG05J 
may be interpreted in terms of nonarchimedean amoebas, by considering the hyper- 
surface defined by the characteristic equation (the eigenvalue A and the entries of the 
matrices A, appearing in the matrix polynomial being thought of as variables) . 

Results of tropical geometry suggest that amoebas obtained by taking the val- 
uation z I— )> log|z| on the field of complex numbers C can be approximated by 
non-archimedean amoebas. Hence, somehow similar results to the ones of [ABG98[ 
IABG041 IABG05| can be expected to hold when using this latter valuation, instead of 
the usual non-archimedean valuation, with equalities replaced by bounds, due to the 
approximation nature. Indeed, the above bounds of Hadamard, Ostrowski, and Polya 
can be seen as the (archimedean) analogue of the results of JABG04J in the particular 
one dimensional case (n = 1). Moreover, in the case of a matrix with non-negative 
coefficients, Friedland j Fri86| established a bound for its spectral radius (or its Perron 
eigenvalue) in terms of a certain maximal circuit mean, which can be interpreted as 
a bound of the maximal eigenvalue of A by the maximal tropical eigenvalue of the 
valuation of A [ABG06,. 

Here we rather replaced the valuation on C by the "valuation" A i— > log || A|| on 



the ring of n x n matrices over C, which leads to (1.5). The idea of using this valua- 
tion was inspired by several works in numerical linear algebra, which suggested that 
the information of the norms is relevant. For instance, Higham and Tisseur [HT03) 
extended to matrix polynomials the bound of Cauchy (related to the special case 
fc = 1 in the Hadamard-Ostrowski-Polya inequality), by using the norms of Aq^Aj 
and A'^^Aj. Fan, Lin and Van Dooren |FLVD04] introduced a scaling based on the 
norms of the matrix coefficients of a matrix quadratic polynomial. In [GS09| . the 



tropical polynomial of (1.5) was initially introduced to refine the results of |FLVD04] . 
It was shown there (see also jShallj ) that the tropical roots of tp can be used to 
perform scalings allowing one to improve the backward stability of eigenvalue compu- 
tations for a matrix quadratic polynomial P {d — 2). Some bounds on the modulus 
of the eigenvalues of P, involving the tropical roots of tp, also appeared in |GS09) 
in the case in which d = 2, and in [Shall] for the smallest and largest tropical roots 
when d ^ 2, which may be seen as particular cases of the lower bounds of the present 
paper. 

Let us finally point out some further related works. Bini used in |Bin96| what we 
call the tropical roots (from the Newton polygon technique) to initialize the Aberth 
method of computation of the roots of a scalar polynomial. Also, Malajovich and 
Zubelli applied Ostrowski's analysis to effective root solving jMZOl] . Finally, Bini, 



Noferini, and Sharify used recently the tropical roots of (1.51 to generalize Pellet's 
theorem [BNS12J. 

The paper is organized as follows. The main results are stated in Section [2] In 
Section |3J we recall the construction of the Newton polygon and give details of the 
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definition of tlie tropical roots. In Section [4| we establish the main lower bounds in 
the case of scalar polynomials, since the proof arguments are more transparent in this 
case. The extension to matrices, and the other results of Section [2j are proved in 
Section [5] In Section [6] we provide examples showing the tightness of the bounds. 

2. Statement of the main results. The inequalities we present here depend 
on the choice of a norm || • || on the space of matrices. We shall in particular consider 
the following assumption: 

(Al) |detA| < II All", for all A € C"^". 
Moreover, the normalized Probenius norm, 






will allow us to generalize Polya's inequality. Therefore, the next assumption will also 
be considered: 
(A2) There exists Q,Q' G C"^", such that detQ = dot Q' = 1 and ||QylQ'||nF ^ 
II All, for all A e C"^". 
We shall see in Section [s] that Assumption |(A2)| implies Assumption |(Al)[ and that 
a number of commonly used norms satisfy Assumption |(A1)| or Assumption |(A2)[ 

Our first theorem provides a lower bound involving the number of nonzero co- 
efficients of a matrix polynomial. This bound is specially useful when the matrix 
polynomial is a "fewnomial" . Recall that the definition of the tropical roots can be 
found in Section [3j see also Figure [2?T| 

Theorem 2.1 (Bounds involving the number of nonzero coefficients). Consider 
a matrix polynomial P with degree d: 

P{X)^Ao + Ai\ + --- + Ad\'^, A, eC"^", O^j^d, 

and let Cit ■ ■ iCnd denote its eigenvalues, arranged by non- decreasing modulus. As- 
sume that II • II is any norm on the space of matrices satisfying \{Al)\ that det Aq ^ 
and let c = L^ n^ . Let ai , . . . , a^ be the tropical roots of the following tropical 
polynomial: 

tp(x) = max (||AJ|a;-') , 

arranged in non- decreasing order. Also let mon P denote the number of nonzero mono- 
mials of P. Then, for all 1 ^ k i^ d, we have 

'^'■■■^"'' ^c{L,r , (2.1) 



(ai ...Qf/c)" 
with 

Lk = -^ ■ (2.2) 

mon_r 

Moreover, when || • || satisfies \{A2)\ the constant L^ can be replaced by the greater one: 

Lk = -^= . (2.3) 

vmon_r 
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Remark 2.2. When n = 1 (thus the matrices Aq, . . . , A^j are scalars), any norm 
is proportional to the normalized Frobenius norm which is nothing but the modulus 
map I • I, satisfies \{A2)\ and for which c = 1. Therefore the tropical roots of tp are 



the same for all norms, the best possible inequality (2.1) is 



iCi---ai 



>L. 



ai 



■Oik 



(2.4) 



and this inequality holds with Lk as in (2.3). 

Remark 2.3. If a norm || • || on the space of matrices satisfies [(Al)| but not |(A2)] 
we can still obtain a bound of the form (2.3) by changing the constant c. Indeed, 
by the equivalence between norms on C"'^", for any norm || • || on C"'^", there exists 
a constant 77 > (which depend on n) such that ||A||nF ^ ^ll^ll for all A G C"^". 
Then, the norm obtained by multiplying |j • || by 77 satisfies |(A2)[ hence (2.3). Since 



the tropical roots of tp and 77 tp are the same, we deduce that (|2.3[) holds for 
with c replaced by c/77". However, if |j • || satisfies [(Al)| but not |(A2)] then 77 > 1, so 
that Inequality (2.3) with c/rf^ may be weaker than Inequality ( |2.2[ ): this is indeed 
the case if and only if 77 > VmonP. The same type of conclusions can be obtained 
for the lower bounds that are stated in the next theorems. 

The following theorem provides a lower bound generalizing the lower bound of 
Polya to the matrix case (since, as said in Remark |2.2[ wh en n = 1, the modulus 
map is a norm satisfying f(A2)l and for which Inequality (2.1 ) reduces to (2.4)). Up to 
the constant c, the following bounds are independent of the coefficients of the matrix 
polynomial P. 

Theorem 2.4 (Universal bound). Lei Aq,..., Ad, P, Ci,---Xnd, \\ ■ \\, c, and 
ai, ... ,ad be as in the first part of Theorem \2.1\ Then, for all 1 ^ k ^ d, Inequal- 



ity (2.1 1 holds with L^ defined as follows: 



8{k){k + l) 



(2.5) 



where £ is defined as in (1.4). Moreover, when \\ ■ \\ satisfies |(A2)[ Inequality (2.1 1 
holds with the greater constant: 



L,= 



v^mk+i) 



(2.6) 



The lower bound of Polya, and its matrix version given above, are tight only for k 
small. By symmetry, one can obtain a tight lower bound when k is close to d (this was 



already noted by Ostrowski in the scalar case |Ost40a| ). Theorem 2.5 below will allow 
us to obtain a tight lower bound when k lies "in the middle" of the interval [0, d], by 
using the comparison between the tropical roots. Unlike the lower bound of Polya, 
this theorem gives a bound which is not anymore independent of the coefficients of 
the polynomial P, although it depends only on a small information, namely a ratio 
measuring the seperation between some tropical roots. Thus, the bound will involve a 
coefficient U{k, 6) depending both on the index k and on the parameter S (the ratio). 
This coefficient is defined as follows: 



U{k,5) :=f(fc) 



1 + 5 + v'A;2(l-5)2 + 4(S 



1-S 



for yfc > and ^ (5 ^ 1 



(2.7) 
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where £{k) is defined by (1.4) for k^ I, and £{0) := 1, and with the convention that 
1/0 = +00, so that U{k, 1) ~ +oo. It is easy to check that 



k+l±l^U{k,S)<e k+^ + ^^ 



1-S 



1- Vs , 



(2.8) 



We shall consider specially situations in which S is small. Then, the following asymp- 
totic regime should be kept in mind: 



U{k,6)r^£(k){k + 1) <e(fc + l), 







Theorem 2.5 (Master lower bound) . LetAo, 



Ad, P, Ci,---,Cnd, II -11, c, and 
ai, . . . , a^ be as in the first part of Theorem 2.1 and denote ao = and oid+i = +oo. 

— '^^ < 1 and 5a. 



Let O^fc <fc<fc"'"^d+l, and denote 5- 



Then, Inequality \2.1\ holds with 



max(L±), L+ :-- 



U{k+ -k-l,5+) 



, L 



with the convention that l/oo = 0. Moreover, when \\ 
Lk of (2.1 1 can be replaced by the greater constant L*^. 



LI. :=max(L 



LV 



1 



LI 



U{k+ 



l,5X) 



^y^ < 1. 



(2.9) 



^ ■ Uik-k-,6-) ' 

satisfies |(A2)[ the constant 



1 



U{k-k-,5l) 



(2.10) 



Theorem 2.5 generalizes Theorem 2.4 and thus the lower bound of Polya. Indeed, 
taking fc~ = 0, fc+ = d + 1, and using that U{k ,0) = £{k){k + 1), we get that the 



constants L^, ^ Lk and Lj 



scalar case {n = 1). 
Note that when (5_ 



^ L^. of Theorem 



2.5 



are exactly the constants Lk of (2.51 



and (2.6) of Theorem 2.4 respectively. Moreover, Theorem 2.5 is already new in the 



1 in Theorem 



2.5 



hl{k — k ,(5_) = +00, so that L^, 







and Lk = L^. Similarly when 5+ = 1, we get L^ = 0. Hence, if (5_ = 1 



5+, we get 
Lk = so that (|2.1|) does not provide any information, although it is true. However, 



if for instance (5_ = 1 and (5+ < 1, we get that Lk = Lt > 0, which gives a positive 



lower bound in (2.1) 



Applying Theorem 2.5 to the particular case when k = k +l = fc+ — 1 lead to 



the following formula for the constants of (2.9) and (2.10) 



and 



max(L±), L+ 



L^=max(Lr), L 



-+* 



1 



Lk* 



l-(5_ 




4(1 + 5_ 


-) ' 


Wi- 


61 



2y 1 + si 



(2.11) 



(2.12) 



But in this case, one can obtain the following stronger lower bounds. 

Theorem 2.6. Let us use the notations of Theorem \2.^ and assume that k = 
fc" + 1 = fc+ — 1. Then, the statements of Theorem 2.5 hold with the constants Lk 



Log-majorization of the moduli of the eigenvalues of a matrix polynomial by tropical roots 7 



and L^, (given in (2.9) and (2.10), or (2.11) and (2.12)j replaced respectively by the 
greater constants L\. and L*^, given by: 



rt) ._ 



1 - 5^64 



(1 



(2.13) 



and 






1 - 6ldl 



(2.14) 



We next indicate how the indices fc+ and k should be chosen, for each 1 ^ k ^ d, 
in order to get the best lower bound Lk- 

Let ko = 0,ki, . . . ,kq = dhe the sequence of abscissa of the vertices of the Newton 



polygon of tp{x), as shown in Figure 2.1 (details on the construction of this polygon 
can be found in Sectionpl. For j = 1, . . . ,q, we have, akj_i+i = • • ■ = ak^ < afc^+i, 
with the convention ad+i = +00. We also denote by 



O^k- 

akj+1 



(2.15) 



for j — 0, . . . ,q, the parameters measuring the separation between the tropical roots, 
in particular 60 = Sq = 0. 



logafc^ 




' log Qd 



fco = ki k2 kj-i kj kq-i kq = d 

Fig. 2.1. Newton polygon corresponding to ip{x). 



Proposition 2.7. Each of the constants L^ and L^, appearing in TheoremH 



maximized by choosing k~ = kj. for some ^ r ^ q such that k^ < k and k'^ — kg + 1 
for some ^ s ^ q such that kg ^ k. This proposition shows that, to apply 
Theorem 2.5 we may always require k~ and A:"*" — 1 to be abscissas of vertices of the 
Newton polygon. Optimizing the choice of k^ , we readily arrive at the following 
corollary. 

Corollary 2.8. Let us use the notations of Theorem \2.5l and let fco = 
0, /ci,...,fcg = d be the sequence of absciss^E of the vertices of the Newton polygon 
oftp{x), as shown in Figure 2.1 Then, the statements of Theorem 2.5 hold with the 



constants Lk, Lf, , L^. and Lj^* (given in (2.9) and (2.10)J replaced respectively by the 
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following optimal ones: 



lT 



- max(i 



±,opt\ 






L--* := 



1 



^y 



o--kj-i<ku{k-kj^i., 



and 



Lr''*:=niax(Lr'°n, 



^k 



(2.16a) 
(2.16b) 

(2.16c) 



(2.17a) 
(2.17b) 



-^^opt ^_ 



(2.17c) 



r.>^.»^^U{k,-k,{^J^) 

1 

_^ . : max . ^^ . 

'■■'^-^<',Juik^k~J^^^) 

However, a simpler choice of fc* consists in taking the nearest vertices of the 



Newton polygon in (2.16) and (2.17), which lead to the following corollary. 
Corollary 2.9. Let us use the notations of Theorem 2.5. let feg = Oj ^ii 



. . , kq 



d be the sequence of absciss <b of the vertices of the Newton polygon oftp{x), as shown 
in Figure 2.1. and let Sq, . . . ,Sq be defined by (2.15). For ^ k ^ d, let us consider 
the unique j G {1, . . . , g} such that fcj-i < k ^ kj, so that ak — a/c _i+i ^ ctk ■ Then, 



the statements of Theorem 2.5 hold with the constants L^, Lj^, L*^ and L^,* (given 
in (2.9 I and (2.10)J replaced respectively by the following ones: 



J prox 

^k '■ — 



f T zli,prox\ 



:t,prox\ 



^k '■- 



U{k,~k,5,) 



J-—, prox 

^k '— 



U{k - /cj_i,(5j_i) 



and 



T+^^prox 

^k '■- 



1 



L*;P™" := max(L±*'P™"), 

1 — ^.prox 

^^' ^k ■ ■ 



1 



/U{k,-k,S]) " ' ^Uik-k,_i,S]_^) 

Moreover, in the particular case where k — kj, we have 

l-^/J7 ' 



(2.18a) 
(2.18b) 

(2.19a) 
(2.19b) 



J- prox \^ i--\-.prox 

^k ^ ^k 



1 



J- ^ ,prox ^ J- +^, prox 
^k ^ ^k 



'1 



1 



Remark 2.10. When all the ratios 8\,...,b„ a re sm all, the maxima in (2.161 



£; 



,opf" 



and (2.17) are attained by taking j as in Corollary 



Li'P™" for all < A: < d. 



2.9 



that is L°^ 



L'^° and 



Remark 2.11. Since U{k,S) is increasing in k and 5, the maximizing j in the 
definition of Lj, *'°^ arises from a compromise between keeping kj ^ k close to k and 
6j small. In particular, when k belongs to an edge of the Newton polygon such that 
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several consecutive edges have almost the same slope than this edge, the maximiz- 
ing j may be the one corresponding to the first vertex at which the slope changes 
significantly, i.e., the first one such that 6j is small. Similar considerations apply to 

7- — *,Opt 

The previous results provide lower bounds for the ratio between minimal products 
of modulus of eigenvalues and minimal products of tropical roots. We next state a 
reverse inequality. 

Theorem 2.12 (Upper bound). Let \\ ■ \\ be any norm on the space of matrices 
C"^". For all i — 1, . . . ,n, we denote by A^"^' the i-th column of A, by rji the least 
positive constant such that ||A*-*'||2 ^ '7i||^|| for all A G c^x" where \\ ■ II2 is the 
Euclidean norm of C", and we set rj :~ rji ■ ■ ■ rjn. Let Aq,. . . , Ad, P , Ci: • • ■ j Qid, c, 



tp, «!,..., ad, be as in Theorem 2.4 let ftQ — 0,ki, . . . ,kq = d be the sequence of all 
absciss(E of the vertices of the Newton polygon of tp{x), and define Sq, . . . ,Sq be as 
in ( |2.15[ ). Denote by Cn.d,k the number of maps : {1, . . . , n} — >■ {0, . . . , d} such that 
Ej=i 4'{j) = k, so that Cn,d,k ^ (min((i, k, nd — k) + 1)""-'^. 



Then, for every j = 1, . . . ,q, if 

I det Afc, 



\A, 



{Cn.d.nk, - l)vSj > , (2.20) 



ICi • • -Cnfejl c f nd 



(ai . . .afc-)" Cj \nk 



Condition (2.201 in this theorem holds if the matrix Ak- is nonsingular and if the 



tropical roots Ukj and afc^^^i are sufficiently separated, so that Sj = akJoi.kj+-i ^ 1. 



Using the results of Corollary |2.9| and Theorem 2.12 we are able to show that the 
modulus of a group of eigenvalues corresponding to an edge of the Newton polygon 
of tp is bounded from above and below by the corresponding tropical roots. The 
following immediate corollary gives an example of such bounds. 



Corollary 2.13. Under the assumptions and notations of Theorem \2.1S\ to- 
gether with (A2) , we have, for all 1 ^ j ^ q, such that Cj > and Cj_i > 0, 



/ nd ^ 

\nkj-iJ 



1-S,Y^' \Cnk,^, + l...Cnk,\ ^ CO /l+J.-lV/' 



.^■■■-,-. •- — ^::^^ <- 



1 + SjJ (afcj"('=^-'^j-i) Cj \l-5j_i 



Similar results can be written when || • || only satisfies |(Al)[ and also for general 
1 ^ k ^ d. The tightness of the bounds in the above corollary depends on the 
parameters Sj-i,Sj and on the condition number of the matrices Ak_i,Ak-. Note 
also that these bounds do not depend on the constant c, except if j = 1 (in which case 
Cj-i = Co = c). Hence, by an argument of continuity, the assumption that det Aq ^ 0, 



which is present there (since the notations and assumptions of Theorem 2.12 include 



the ones of the first part of Theorem 2.1 ), can be dispensed with, except when j = 1. 



3. Tropical polynomials and numerical Newton polygons. We recall here 
basic results on tropical polynomials of one variable. See for instance |BCQQ92llVir01[ 
IIMS07| for more background on tropical polynomials from different perspectives. 

Let Mmax denotes the set M U {— c»}. A (max-plus) tropical polynomial / is a 
function of a variable x € Umax of the form 

f{x) = max (/j +jx) , (3.1) 
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where d is an integer, and /o, . . . , /d are given elements of M„iax. We say that / is of 
degree d \i f^ ^ — oo. We shall assume that at least one of the coefficients fo, ■ ■ ■ , fd 
is finite (i.e., that / is not the tropical "zero polynomial"). Then, / is a real valued 
convex function, piecewise affine, with integer slopes. 

Cuninghame-Green and Meijer showed |CGM80] that the analogue of the fun- 
damental theorem of algebra holds in the tropical setting, i.e., f{x) can be written 
uniquely as 

d 

f{x) = fd + '^ma,x{x,aj) , 

j=i 

where cti ^ • • • ^ a^ G I^max- The numbers ai, . . . ,ad are called the tropical roots. 
The finite tropical roots can be checked to be the points at which the maximum in 
the expression ( |3.1[ ) of f{x) is attained at least twice, whereas — oo arises as a tropical 
root if /o = — cxD. The multiplicity of a root a is defined as the cardinality of the set 

{j e {l,---,i} I aj = «}■ 

The multiplicity of finite root a can be checked to coincide with the variation 
of the derivative of the map / at point a, whereas the multiplicity of the root — cxd 
is given by inf{j | fj ^ —00} or by the slope of the map / at —00. The notion of 
tropical roots is an elementary special case of the notion of tropical variety which has 
arisen recently in tropical geometry |IMS07| . 

The tropical roots can be computed by the following variant of the classical New- 
ton polygon construction. Define the Newton polygon A(/) of / to be the upper 
boundary of the convex hull of the region 

{{j, A) e N X M I A s^ /j, s^ j ^ rf} C K^ . 

This boundary consists of (linear) segments. 

The following result was established in |ABG05| . it relies on standard Legendre- 
Fenchel duality argument. 

Proposition 3.1 ^ |ABG05[ Proposition 2.10]). There is a bijection between the 
set of finite tropical roots of f and the set of segments of the Newton polygon A(/); the 
tropical root corresponding to a segment is the opposite of its slope, and the multiplicity 
of this root is the length of this segment (measured by the difference of the abscissiE of 
its endpoints). (Actually, min-plus polynomials are considered in |ABG05] . but the 
max-plus case reduces to the min-plus case by an obvious change of variable.) 

Since the Graham scan algorithm |Gra72| allows one to compute the convex hull 
of a finite set of points by making 0{n) arithmetical operations and comparisons, 
provided that the given set of points is already arranged by abscissae, it follows that 
the tropical roots, counted with multiplicities, can be computed in linear time (see 
also |GS09[ Proposition 1]). In particular, the maximal tropical root is given by 



ad — max 



fj - Id 
CKd d~ j 

Example 3.2. Consider 

f[x) = max(0, 1 + .T, 6 + 2a;, 4 + 4a;, 9 -I- 8a;, 5 + 10a;, 1 + 16a;) 



The graph of / and the Newton polygon of / are shown in Figure |3.1| The tropical 
roots are —3, —1/2, and 1, with respective multiplicities 2, 6 and 8. 
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Fig. 3.1. Graph of the tropical polynomial function f on left, and of its Newton polygon on right. 



The notion of root also applies with trivial changes to the "max-times" model of 
the tropical structure, in which polynomial functions now have the form 



tp{x) 



max ajX-' 



where aQ,...^a^ are nonnegative numbers, and the variable x now takes nonneg- 
ative values. Then, the tropical roots of tp{x) are, by definition, the exponen- 
tials of the tropical roots of its log-exp transformation f{x) :— logtp(exp(a:)) = 
maxo^j^d(logaj- + jx). 

In the sequel we shall consider max-times polynomials associated to usual scalar 
or matrix polynomials. We shall need the following result which follows from the 
above definitions and properties. 

Proposition 3.3. Let tp{x) = maxo^j^d Q-j^;-' be a max-times polynomial of 
degree d, with Oj ^ 0, j = 0, . . . d. Assume that oq 7^ 0. Let < ai ^ • • • ^ a^ denote 
the tropical roots oftp arranged in non- decreasing order. Then, 



aj < ooj^a^ \ 



0, 



e=i 



Moreover, let k he the abscissa of a vertex of the Newton polygon of p, then Ok > 
and 

aj < flfc nLfe+i "7^ forj^k,...d , 
aj < flfe rifcj+i "^ forj = 0,...k . 

Proof. Let f{x) := maxo^j^d(logaj + jx). By definition logofj, j — 1, . . .d, are 
the tropical roots of /. The upper boundary of the Newton polygon of / coincides 
with the graph of the concave hull a of the map j G {0, . . . ,d} 1— > logOj e M. By 
Proposition |3.H the tropical roots are the opposites of the slopes of d and their 
multiplicities are the lengths of the segments where d has this slope. This means that 



log aj 



ij, hence Oj 



a-o 



X^fci logof^j ^i^d using that d is above the map 



J 1-^ logaj, and that both maps coincide at the boundary point j = 0, we get the 
first inequality of the proposition. If now k is the abscissa of a vertex of the Newton 
polygon of p, then {k,ak) is an exposed point of the Newton polygon, which implies 



that dfc = logflfc. Since dj = d^ 
inequalities of the proposition. D 



Sfcfc+i logct^i foi' ^-U J ^ ^; ■^^ get the two last 
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4. Tropical bounds for the modulus of the roots of a scalar polynomial. 

In this section, we prove the main lower bounds of Section [2] for scalar polynomials, 
since the arguments are more transparent in this case. The generalization to the 
matrix case, as well as the proof of the upper bound, will be given in the next section. 



From Remark 2.2 in the scalar case, all lower bounds reduce to (2.4) with the 



constant Lk obtained under Assumption (A2) that is (2.3) in Theorem 2.1 (2.6) in 



Theorem 2.4 (2.10) in Theorem 2.5 and (2.14) in Theorem 2.6 The proof of these 



lower bounds is based on the next result proved by Landau in |Lan05] . building on 
an earlier observation of Lindelof |Lin02j . This result is one key step of the proof of 
Polya's inequality [Ost40a]. 

Lemma4.1 f [LanOSp . LetC,i,...,C,d he the roots of a polynomial p{x) ='^j^qO,jX^ , 
arranged by non- decreasing modulus, and assume that oq 7^ 0. For all I ^ k ^ d, we 
have 



logic^^ai^-infiiogftjai; 



■J I ^2a-fe) 



(4.1) 



We include the short proof, as its idea will be used in the extension to the matrix 
case given in the next section. 

Proof. The formula of Jensen |Jen99j shows that if Ci, . • . , Cfe ^re the roots of p(z) 
in the closed disk of C of radius r, counted with multiplicities, then 

log ICi . . . ai = fclogr + log b(0)| - ;^ / " log \p{re'')\d9 . 
It follows that, for all r > 0, and k = 1, . . . ,d: 



log|Ci...ai ^A:logr + logb(0)| 



27r 



log|p(re*^)|d6' 



(4.2) 



Using the comparison between the geometric and the L^ mean, together with Parse- 
val's identity, we get 

^ fj log \p{re^')\de ^ \ log {^ J^^ \pire^')\'de) 



- log ( — 

2 ^^TT 



2Tr d 



\J2a,r=e^'')\'dO) 



j=o 



1 " 

ilog(5:|a,|V^) 



j=o 



Gathering this inequality with (4.2) yields 



1 " 

log|Ci...ai^fclogr + log|ao|--log(^|a,|V^) 



3=0 



Since this holds for all r > 0, this shows the inequality of the lemma. D 

Using the same arguments as in the proof of Polya's inequality reproduced in |Ost40a| , 
we obtain the following result which involves now the constants ak instead of the mod- 
ulus jofcl of the coefhcients of p. 
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Corollary 4.2. Let ^i, . . . , C^ be the roots of a univariate scalar polynomial p 
of degree d, p{x) = X]i=o S^;-', arranged by non- decreasing modulus, and assume that 
ao 7^ 0. Let ai, . . . , a^ denote the tropical roots of the associated tropical polynomial 
tp{x) = maxo^j^d |aj|a;^, arranged in non- decreasing order. Then, for all 1 ^ k ^ d, 
Inequality (2.4) holds with L^ such that: 



{LkY 






?> 



(4.3) 



with 



Pk, 



ifi = k 



(4.4) 



which satisfies in particular (3k,j ^ 1 for all k,j — 0,...,d. 

Proof. By applying the change ol variable r = ctkV^ in the inequality ol Lemma 
we get 

iogict...ai.sup(4iog(EM(.k)-)) 



4.1 



«>o 



,£=1 



\J=0 



Applying Proposition 
for all j = 0, . . . , d, w 



3.3 



to the max-times polynomial tp, we get \aj\ ^ |ao| Ilfci '^e 



rich with the above inequality yields (2.4) with 



{LkY 




iuli^^mi 



^')^ 



k-j 



CJ-k 



which can be written in the form (4.3) with (S^.j as in (4.4). Moreover, since a, is 



nondecreasing with respect to j, we get that all the Pk.j sue less than or equal to 1. D 
Polya's inequality is then an immediate consequence of Corollary |4.2[ Indeed, 
using the property that all the /3k, j are less than or equal to 1, we obtain: 



{LkY 



< 



d 
J=0 



:j-k 



^E^^-'- 



1 



e(i-0 



for all ^ > 0. The minimum of the right hand side of the previous inequality for 
< ^ < 1, is attained for ^ — k/{k + 1), from which we deduce (1.31, which is also 
the scalar version of Theorem 12.41 



We can now deduce similarly, from Corollary |4.2[ the scalar versions of the lower 
bounds of Theorems |2.1[ |2.5[ and |2.6[ Since in the scalar case, we are reduced to 
show (2.4) for some constants Lj., and that this inequality is precisely the statement 



of Corollary |4.2[ we only need to show that these constants Lk are lower bounds of 



the constant Lk of (4.3) 
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Proof of the scalar version of Theorem 2.1 Using the property that /3fej ^ 1 for 



all fc,j, we obtain that the constant Lk of (4.3) satisfies {Lk) ^ ^ X) 



i=0,aj#0 



all f > 0. When ^ = 1, the right hand side of this inequality is equal to the number of 



non-zero coefficients of p, which shows that the constant Lk of (4.3) is lower bounded 



by the constant Lk of (2.3), which implies (2.4) with this lower bound Lk. D 

Proof of the scalar version of Theorem 2.5 Assume that ^ fc^ < /c < fc+ ^ d+1. 
Denote 5- = ak-joik and (5+ = ak/cik+- Then, for all 1 ^ £ ^ fc, a^ ^ a^, and for 

all 1 ^ £ ^ fc^, Q!£ ^ ak- — ctk^-- This implies that Pk.j ^ 5L "', for all j ^ k^ . 
Similarly, for all i ^ k, ai ^ ak, and for all £ ^ k^ , ai ^ ak+ = ak/S^, hence, for all 

. Since we also have /3k,j ^ 1 for all j, fc, we obtain that the 



j^k+, pk,j < <5i 



constant Lk of (4.3) satisfies, for all 5^ < ^ < 5, ^, ^ 7^ 1, 



fc+-i 



yJ=0 / \j=fe-+l 



E^' 



2(]-k++l)gj^k 



=k+ 



<? 



k^-k+1 



1 tk -fe+1 _ fA;+-fc 



-k-1 su 



^~5l 1-e 

Since "t2 ^ ^ ^^2 c ^he later inequality can be written 

(Lfc)-2<g(0, V<5!<^<<5+2, ^^1 



1-ne 



as 



(4.5) 



with 



9{i) 
9-iO 



g~{0+9+i0 , 



-fe+1 



1 



5+(e):=? 



fe+-fe-i 



1 



1 



Note that (4.5 1 also holds when fc~ = or /s+ = d + 1, since then (5_ = or 5+ = 
respectively. When 5- = 1 and J+ < 1, the conditions on ^ in (4.5) are equivalent to 
1 < ^ < (5^^, whereas when (5_ = (5+ = 1, these conditions are never satisfied, but in 



this case the constant Lk of (2.10) is equal to so there is nothing to prove. 

The functions g_ and 3+ satisfy 5_(C) = 9k~k- ,sl{i~^) andg+(0 = gk+~k-i,sl{0 



where gk,5 is defined for ^ 7^ 1 and 5 ^ by: 



_ ffe+i 



1-6 



We have gk.siC) < for ^ < 1, hence 5_(0 s^ for aU ^ > 1, and 5+(0 < for all 
^ < 1. When (5 < 1, the minimum of gk,5 on (1, 5^^) is attained at 



ik. 



s ■- 



k(S+l)-^k'^(S+lY-iS(k^-l) 
2(fe-l)<5 

2 
1+5 
fe + 1 

fe 



when fc 7^ 1, 5 7^ , 
when fc = 1 , 
when 5 = 0. 



The last formula gives £^k.s = 00 when (5 = and fc = 0, which is the point of infimum 
of 5fe.5 for 1 < ^ < 00 = (5^^, since gk,s is decreasing. It also gives ^^.i ~ 1 for all fc ^ 0. 
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Hence extending gu.s, .9- and g^ by +00 at point 1, we get that the infimum of g^^s 
on (1,(5^1) equals gkMik,s), and denoting ^_ = Hk-k-j^y^ and (,+ = ^k+~k^i,sl, 
we obtain 



(Lfe)-2 ^ min(g(^_),5(C+)) s^ mm{g_{^_),g+{^+)) . 

In order to simplify this bound, we need to find good estimates of £,k,s- 
When k ^ 1 and J 7^ 0, we have: 



(4.6) 



^k,S 



k{6+l)-y/k^il-6y+4S 
2{k-l)5 

2(fc + l) 

' k{6 + 1) + ^fc2(l-5)2+4j ■ 

Using the property that ^ ^ in the last formula, we get that £,k,s ^ ^^ = £,k,o- 
Moreover, this inequality also holds for k = 1. In particular {^k,s)'' ^ £{k), for all 
k ^ and S ^ (taking the convention that ^" = 1 for all ^ e (1, 00]). We also have, 
for fc ^ and ^ (5 < 1, 



1 



This yields 



1 



1 - 5^k,s 



l + S+Jk^(l-S)^+AS , 1 

: ^ k ' 



V~5 



1-Vs 



gkA^k,s)^U{k,S) Vfc^O, 0^,5s^l, 



with U as in (2.7 1. Indeed, this inequality holds for (5 < 1 by the above arguments. 
It also holds for 6=1 since U{k, 1) = +00. This implies in particular that 5_|_(^+) ^ 
I4{k~^ — k — 1,(5^) and g_(^_) < U{k — k^,S^). Combining these two inequalities 



with ( |4.6[), w e obtain that the constant Lk of (4.3) is lower bounded by the constant 
L^. of (2.10), which implies ( |2.4[ ) with this lower bound L^. instead of Lfc. D 

Remark 4.3. In the previous proof, obtaining the minimum of g instead of g_ 



and g+ would have led to a better lower bound for Lk than the one of Theorem 2.5 
However, such a bound is more difficult to estimate for general values of k^ — k~ . We 



can use for instance the first inequality in (4.6), which will give better estimates in 
some particular cases. For instance, when fc = 0, we get that gk,s{£,k,s) = U{k,8) = 
(1 + \/(5)/(l — v^). In particular, when fc = fc+ — 1, the bound {Lk)~^ ^ g+{^+) gives 
the lower bound L'^ of (2.12), with 5 = 5^. However (4.6) gives the slightly better 
bound: 



(Lk)-^ ^ g{^+) = 



1 + S+ 1 - (52 



1-6+ (l - 5+5^_)il - 5+) 



Proof of the scalar version of Theorera \2.6\ Let us use the notation of the previous 
proof. When fc = A:+ — l = fc^ + l, Inequality ( |4.5[ ) is true for all ^ G {5'^,5+'^) since 
g can be well defined at 1 (by continuity for instance). When 5^5+ < 1, the interval 
{S^ ,S+ ) is nonempty, the map g is convex there and its minimum is achieved at the 



point^ = {l + 6lS+)/{6+{l + 5+)). Then ([45]) yields (Lfc)-2 ^ (1 + 5+)'^ /{I- 6^1) 
which implies that the constant Lk of (4.3) is lower bounded by the constant LV 



of (2.14), which implies (2.4) with this lower bound L?, instead of Lk- Moreover, 
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since the minimuni of g is les s or equal to the right hand side of (4.6), the proof of 
the scalar version of Theorem 



2.5 



implies that L*/^ ^ L^.. 



U 



Example 4.4. Consider the following scalar polynomial 

p{x) = l-exp(l)a;-exp(6)a:^ + exp(4)x'* + exp(9)a;^+exp(5)a;^" + exp(l)a::^'^ . (4.7) 

The log-exp transformation of its tropical polynomial tp is the tropical polynomial / 
The graph of / and the associated Newton polygon were shown in 
Hence, the tropical roots of tp are the exponentials of —3, —1/2, and 1, 



of Example 3.2 



Figure 3.1 



with multiplicities 2, 6 and 8, respectively. In Table |4T| we are comparing the value of 
the ratio ^''^'"''''^ , with the lower bound of Polya given in (1.3), together with the lower 



bounds obtained in Theorem|2.1|and Corollary|2.9| There, for each k = 1 



the column RATIO gives 



ICl 



16, 



the column POLYA gives the lower bound of Polya, 



ai...ak 

which coincides with the constant Lk of (|2.6[), the column fewnom gives the lower 



bound Lk of (2.3), and the column SEPAR gives the lower bound 7^*;?''™ of (2.19) 



based on the separation between the tropical roots. 



k 


RATIO 


FEWNOM 


POLYA 


SEPAR 


1 


0.93475 


0.37796 


0.5 


0.5 


2* 


1.00034 






0.3849 


0.92102 


3 


0.98782 






0.32476 


0.49664 


4 


0.97926 






0.28622 


0.38360 


5 


0.98348 






0.25880 


0.32403 


6 


0.98771 






0.23802 


0.37508 


7 


0.99383 






0.22156 


0.47570 


8* 


1.00000 






0.20810 


0.79696 


9 


0.98811 






0.19683 


0.47570 


10 


0.97636 






0.18721 


0.37508 


11 


0.96476 






0.17888 


0.31917 


12 


0.95329 






0.17158 


0.28622 


13 


0.96476 






0.16510 


0.32476 


14 


0.97636 






0.15930 


0.3849 


15 


0.98811 






0.15407 


0.5 


16* 


1 






0.14933 


1 



Table 4.1 

Comparison of the lower bounds for the scalar polynomial p of 
vertices of the Newton polygon oftp are indicated by the symbol *. 



4.7[l. The abscissa; of the 



Example 4.5. Consider now the following scalar polynomial 

p{x) =1 — exp(3)a; — exp(6)a;^ — exp(7)a;^ — exp(8)x^ + exp(9)a;* 
+ exp(7)a;^° - exp(4)a;" - exp(3)a;^'' + exp(l)a;^^ . 



(4.8) 



The log-exp transformation of its tropical polynomial tp is again the tropical poly- 
nomial / of Example 3.2 but now there are some points of the graph of the map 
k I— >■ log \pk I , where pk is the fcth coefficient of p, that are on the edges of its New- 
ton polygon, that is its concave hull. Table |4.2| is obtained on the same principle 
as Table |4.1| In Figure |4.5[ we plot in the same graph the values of the ratios and 
lower bounds for the scalar polynomials p of (4.7) and (4.8). This indicates that the 
lower bound Lk of Corollary |2.9[ based on the separation between tropical roots, may 
become tighter when the polynomial has non zero coefficients between the vertices of 
the Newton polygon. 

5. Proofs of the bounds for the modulus of the eigenvalues of a matrix 
polynomial. In this section we provide the proofs of the main results which were 
stated in Section [2] and check the properties given at the beginning of Section [2] on 
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k 


RATIO 


FEWNOM 


POLYA 


SEPAR 


1 


0.61759 


0.31623 


0.5 


0.5 


2* 


0.98684 






0.3849 


0.92102 


3 


0.84948 






0.32476 


0.49664 


4 


0.73124 






0.28622 


0.38360 


5 


0.63280 






0.25880 


0.32403 


6 


0.54760 






0.23802 


0.37508 


7 


0.72118 






0.22156 


0.47570 


8* 


0.95684 






0.20810 


0.79696 


9 


0.76959 






0.19683 


0.47570 


10 


0.61899 






0.18721 


0.37508 


11 


0.57724 






0.17888 


0.31917 


12 


0.53832 






0.17158 


0.28622 


13 


0.62117 






0.16510 


0.32476 


14 


0.71677 






0.15930 


0.3849 


15 


0.84662 






0.15407 


0.5 


16* 


1 






0.14933 


1 



Table 4.2 

Comparison of the lower bounds for the scalar polynomial p 
vertices of the Newton polygon oftp are indicated by the symbol *. 



of R 



The abscisses of the 




Fig. 4.1. Plot o/ RATIO as a function of k for the scalar polynomials of |4.7[ l (blue circles) 
and ( |4.8[ l (red squares), together with the lower bound FEWNOM (blue and red lines resp,), POLYA 
(plus signs) and separ (stars signs). 



In particular as in Theorem 2.12 
of matrices C 



norms on the space of matrices. For all p G [1, cxd], and n ^ 1, we shall denote by || • ||p 
the P norm of C": \\v\\p := (X]"=i kil'')"'"^'' for p < cx) and ||f ||oo '■— niaxi=i^...^„ \vi\. 

j 2 is the Euclidean norm. The norm on the space 
induced by the norm || • ||p on C" (the same norm is used for the 
domain and the range of A) will also be denoted by || • ||p. In particular the norm 
II • II2 on C"^" is the spectral norm. Moreover, for all p G [l,oo), and n ^ 1, we shall 
denote by || • ||*p the following normalized Schatten p-norm on the space of matrices 
C"^": II A||*p :— (^ X]"=i ^lY^^^ where cti, . . . , cr„ are the singular values of A (the 
eigenvalues of y/A*A, where A* denotes the conjugate transpose of A). Then, the 
normalized Schatten 2-norm coincides with the normalized Frobenius norm. Recall 
that the (unnormalized) Schatten 1-norm is also called the trace norm or the Ky Fan 
n-norm. We shall also denote by ||A||*co '■— nia'Xi=i_....„ cr^ the Schatten cx)-norm, 
which coincides with the spectral norm. 

Property 5.1. The normalized Frobenius norm satisfies Assumption |(A1)] 
Therefore, Assumption\{A'l)\ implies Assumption\{A'C)\ 
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Proof. From Hadamard's inequality, we get that |detA| s^ ||yl(^)||2 • • • ||A("'||2, 
for all A £ C"^", where A^^^ denotes the j-th column of A. Then, since the geometric 
mean is less than or equal to the arithmetic mean (or simply by the concavity of 
the logarithm), we get that |det^p/" < ^(||A(i)||^ + ••• + P^"^Hi) = \\A \\lp, which 
implies that the normalized Frobenius norm satisfies Assumption |(Al)[ Since for 



all Q,Q' e C"^" such that detg = detQ^ 1, we get |dety4| =_dct(g^g')| ^ 
\Q^Q' WnF for all A e C"^", Assumption |(A2)| implies Assumption |(A1)[ D 



s« 



Property 5.2. For all p e [l,oo], the norm \\ ■ \\p on C"^" satisfies (A2) 
Proof. Indeed, let us denote by ||A||p,(j the norm induced by the norms | 
and II • llg on the range and domain C" of A respectively, which means that || A||p_, 
max{||Aw||5/||t;||p, v e C",v ^ 0}. Since (for aU v e C") p G [l,cx)] ^ \\v\\p eR+ 
is a nonincreasing map, we get that HAHp^g ^ ll^llp'.g' when p ^ p' and q ^ q'. 
In particular ||A||p = |l^||p,p > ||^||2,oo when 2 < p, and ||^||p = ||A||p,p ^ ||^||i,2 
when p ^ 2. In addition, it is easy to show that ||A||2,oo = niax{||(A*)(*'||2; J = 
l,...,n}. Since the last expression is greater or equal to ||^||nF, we deduce that 
ll^ljp > PII2.00 ^ PlInF for aU A e C"^" and 2 < p ^ 00. Similarly, |lA|li.2 = 
max{||AW||2,i = l,...n} ^ PlInF, hence ||A||p ^ ||A||i^2 ^ P||nF for all A e C"^" 



and 1 ^ p ^ 2. This shows that || • ||p satisfies (A2) with Q = Q' = I the identity 
matrix. D 

This result implies that all norms induced by the norm ||w||d ~ ||gu||p in the 
domain of A and the norm ||u||r — IIQ'i'lIp in the range of A with the same p, but 
possibly different matrices Q, Q' such that det Q = det Q' = 1, satisfy [("A2)[ 

Property 5.3. // 1| • || is the norm on C"^" induced by a norm on C" (the same 
for the domain and the range of matrices), then it satisfies \{Al)\ 

Proof. Since det A is the product of the eigenvalues of A counted with multi- 
plicities, we get that | det A| ^ p(A)", where p{A) is the spectral radius of A. Since 
p{A) ^ ||A|| for all A e C"^" and all induced norms || • || on the space of matrices, we 
get the result. D 



Property 5.4. The normalized Schatten p-norm \\ ■ ||*p on C"^" satisfies (Al) 
for all p G [1, 00] . It satisfies |(A2)| if and only if p ^ 2. Moreover, for p £ [1, 2), the 



least constant rj such that t]\\ ■ ||h,p satisfies (A2) is given by rj = n^'P ^'^ > 1 



Proof. Since | det A\ is the product of the singular values of A counted with mul- 
tiplicities, we get that | det A\'p = crf ■ ■ ■ a^, and using that the geometric mean is less 
than or equal to the arithmetic mean, we obtain that | det A\p/'^ ^ - (a^ + • • • + cr^) 



yl||^p, which implies that the normalized Schatten p-norm satisfies Assumption (Al) 
We have that p i— > ||^||*p is a nondecreasing map. Hence, when p G [2,cxd] 
A||*p ^ 11^11*2 = II^IInF, for ah A e C"^", which implies that || • ||*p satisfies 
when p ^ 2. Let 77 > be the least constant such that ?7|| • ||*p satisfies 



(A2) 
(A2y 



(A2)| for p < 2, which will finishes 



(A2)| if and only if p ^ 2" . Since 



and let us show that rj — n^/p~^/'^ when p < 2. This will implies in particu- 
lar that 7/ > 1 hence || • ||*p does not satisfy 
the proof of the equivalence "|| • ||*p satisfies 

"-^11 ' II *p is the unnormalized Schatten norm, which is nonincreasing with respect 
to p, we get that, for all p e [1,2), and A S C"^", n^/^pn^p ^ n'^/P\\A\Up, 
which implies that 77 < n^/P-'^^'^. Let us fix p, Q,Q' e C"^" such that detQ = 
detg' = 1 and ||gAg'||nF ^ v\\M*p. for all A e C"^". For alii = l,...,n, let 
us consider the matrix A whose entries are all zero but the entry ii equal to 1. 
Then, the singular values of A are all zero except one which is equal to 1, so that 
Pll^p = {l/n^/P. Moreover g^ig' = g^'^g^^), where Q',^. denotes the i-th raw 

of g'. Hence ||gAg'||„F - {l/n)'/^Q^'^h\\{Q'*)'^'M\2- From ||gAg'||nF < v\\M*p, 
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we deduce that ||Q('^||2||(Q'*)**'||2 s^ ^{\|nfl^-^l'^ . Since detQ = detQ' = 1, we 
get using Hadamard's inequality, 1 s^ ||Q'''||2||(Q'*)(')||2 • ■ ■ ||Q'"^||2||(Q'*)("^||2 < 
(77(l/n)i/P-i/2)n^ ^j^i(,jj gjjQ^g ^j^g^^- ,^ ^ ^iIp-\Ii_ ^^ dcducc that T] = n^/p-^/'^ > 1 

for p > 2, which completes the proof. D 

For the proof of lower bound results stated in Section [5] the approach is similar 
to the scalar case, except that we change the application of Jensen formula to obtain 
the following matrix version of Lemma |4.1| 

Lemma 5.5. Let Ci, . . . , (nd denote the eigenvalues of the matrix polynomial P 
with degree d: 



PiX) = Ao + AiX + ---+AaX'^ A^ e 

arranged by non- decreasing modu lus. Assume that 
of matrices which satisfies 
I ^ k ^ d, we have 



(Al) that det^o 7^ and let c ~ Krj- 



is any norm on the space 
- ^i^^^^ Then, for 



log ICi . . . Cnk\ > logc - n inf log V 



\A 



r>0 



VJ=0 



I An 



lA^]-k 



(5.1) 



When II • II s atisfies\{ A2)\ the previous bound can be improved as follows 



log|Cr...C.|^logc-.infilog(^gM 



.2(i-fe) 



(5.2) 



Proof. From Inequality (4.2 1 applied to p{X) — detP(A), we get, for all r > 0, 



log ICi • • ■ Cnk\ > nfclogr + log I det Aol ~^ \ log | det P(re^'')|d0 . (5.3) 

Since || • || satisfies [(Al)[ we have 

d d 

|detP(re*«)| ^ \\P{re^'W = || ^^ A,(re*«)^|r < ^\\A,\yr , (5-4) 



j=o 



i=o 



for all Q G [0, 27r]. Gathering ( [KI] ) with (|53| , we obtain ([53 1. 

Assume now that || • || satisfies [(A2) with some matrices Q, Q' such that det Q = 
detQ' = 1, a nd let us prove (5.2). Since from Property 5.1 the normalized Frobenius 
norm satisfies |(Al)[ we get that | det P(re'^)| = | det(QP(re''')Q')l < IIQ^(^e*^)'3'linF- 
Now using the comparison between geometric and L? means, we deduce 



^- I log\detP{re'')\d9^n^ f log ||gP(re^'')Q'||n 
'■"' Jo ^^ Jo 



^d9 



27r 



^fiog(^y^ \\QP{r^')Q'\\l^de 



(5.5) 



From the formula of the normalized Frobenius norm, we get by applying Parseval's 
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identity to each coordinate {QP{re^^)Q')im 

£,m=l \j=0 J 

d 

= Y,{\\QA,Q'\\lyr^') 
j=o 



'de 



3=0 



Gathering this inequality with Inequalities (5.5 1 and (5.3), we obtain (5.2). D 



Using the same arguments as for Corollary |4.2[ we obtain the following result. 
Corollary 5.6. Let Ci, . . . , C,nd denote the eigenvalues of the matrix polynomial 
P with degree d: 



P{X) = Ao + AiX 



AaX" A^e 



< J s^ d , 



arranged by non- decrea sing m odulus. Assume that || • || is any norm on the space of 



matrices which satisfies 



(Al) 



that det Aq ^ and let c 



_ IdctAol 



IIAn 



^. Let ai, . . . ,ad be 



the tropical roots of the following tropical polynomial: 

tp(x) ~ max (WAAlx-') 



arranged in non- decreasing order. Then, for all 1 ^ k ^ d, Inequality (2.1 ) holds with 
Lk such that: 



(ifc)~' = inf E P'^^^^- 



j'k 



(5.6) 



,j=0,Aj#0 



with Pk.j as in (4.4). 

Moreover, when \\ ■ \\ satisfies \{A2)\ the constant Lk in (2.1) can be replaced by 
the greater one: 



(ur- 






(5.7) 



Proof of Theorems 2.1 2.4 2.5, and 2.6 In Corollary 5.6 the formula (5.7) oi Lk 



coincides with (4.3), except that the condition Oj ^ is replaced by Aj ^ 0. Hence, 



the lower bounds of the constant Lk of (4.3) obtained in the proofs of the scalar 



versions of Theorems |2.1[ |2.5[ and |2.6| give similar lower bounds for the constant 



Lk of (5.7), which combined with the assertion that Inequality (2.1) holds with this 



constant, now provide the lower bounds of Theorems |2.1[ |2.5[ and |2.6| respectively. 



in the case where || • || satisfies |(A2)[ that is Inequality (2.1) with Lk replaced by 
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and the constant L*j^ of (2.14) 



the constant Lk of (2.3), the constant L*/, of ( |2.10[ ) 
respectively. 

In the case of a norm satisfying only |(Al")| the formula (5.7) oi L^ is replaced 
by (5.6), which means that everything behave as if Lk were replaced by its square 



and the numbers a^ were replaced by their square roots in (5.7). Hence the assertions 
of Theorems 2.1[ 2.5 and 2.6 in the case of Assumption (A2) are still true under 



Assumption (Al) up to this transformation of the constant L^ of (2.3), the constant 



L*f. of (2.10), and the constant L^ of (2.14), respectively. 



The constant Lk of (2.3 ) does not depend on the values of the constants ak, hence 
the above transformation leads to the square of Lk , which is exactly the constant Lk 

The constant LI of (12. 101) is a 



of (2.2), which finishes the proof of Theorem 



2.1 



,d+l, 



function of S^ — -^- and (5+ — -2^^. Since replacing all numbers a,, i = 0, 

by their square roots, reduces to replace the numbers 5^ and 5^ by their square roots 



too, the above transformation of the constant L^ of (2.10) consists in taking its square 



and replacing 5_ and (5+ by their square roots, which leads to the constant Lk of (2.9) 



-*tt 



which finishes the the proof of Theorem |2.5| Similarly, the above transformation of 



the constant L^ of (2.14) leads to the constant L\. of (2.13), which finishes the proof 



2.5 



of Theorem HJl 

Now Theorem |2.4| is an immediate corollary of T heorem |2.5 
Section^ taking /c^ = 0, fc+ = d + 1 in Theorem 
andW(fc,0) =f(fc)(fc+l), we get Theorem [zil D 

Corollaries |2.8| and |2.9| are straightforward consequences of Theorem |2.5| and 
Proposition |2.7| We only need to prove the latter proposition. 



Indeed as said in 
and using L^ < ife, i^* < L^. 



Proof of Proposition 2.1 Let us use the notations of Theorem 2.5 and let fco = 
0, /ci, . . . , fcg = d be the sequence of abscissas of the vertices of the Newton polygon of 
\.p{x), as shown in Figure 2.1 Then, 
inside an interval [fcr-i 



maximize L, or L, 



the ratio (5_ does not change when k~ moves 
1, kr\. Since U{-, 5) is a nondecreasing function, in order to 
and thus i^, with 5- constant (and 5+ and fc"*" constant), we 
need to minimize k ~ k^ , which implies that k~ = kr for some ^ r ^ g such that 
kr < k. Similarly, fc+ = fc^ + 1 for some ^ s ^ q such that kg ^ k. D 

Remark 5.7. AU the lower bounds in Theorems |2.H |2.4[ |2.5[ and |2.6| are equiv- 
alent to inequalities of the form 



(|Ci...Cnfcr'|detAo|)'/"ifeSSPo||ar' 



.a,. 



(5.; 



By definition of the tropical roots, || Ao||aj~ . . . a^T = expp^, where pk is the value in 
k of the concave hull of the map j G {0, . . . , d} i— )• log \\Aj\\. Hence, if |j • || and || • |j' are 
two norms on the space of matrices such that ||^|| ^ ||^||' for all A g C"X"^ then the 
right hand side in ( |5^ is smaller for || • || than for || • ||'. Moreover, in ( p^ , ( [231 ), PiU.' 
and (2.6), Lk depends only on P or fc and n, but not on the norm || • |j, so that ( |5.8[ ) 
is necessarily a tighter inequality for || • || than for || • ||'. In particular if |(A2)| holds 
with Q and Q' , then the lower bounds (2.3) and (2.6) for || • || are weaker than the 
corresponding ones for A H> ||QAQ'||nF. For the lower bounds (2.9), ( 2.10[ ), ( |2.13 ), 
and (2.14) of Theorems 2.5 and 2.6 and the ones of Corollaries 2.8 and 2.9 based on 



the separation between the tropical roots, the comparison is not so simple, because 
changing the norm changes the separation between the tropical roots. 

To prove Theorem |2.12[ we first prove the following lemma which provides a lower 
bound for the modulus of the coefficients of the polynomial det(P(A)). 

Lemma 5.8. Let Aq,. . . , Ad, P, \\ ■ \\, rj, ai, . . . ,ad, q, 5q, . . . ,6q, and Cn,d,k be 
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as in Theorem \2.1S\ and denote by p the polynomial: 



nd 

p(A)=detP(A) = ^p,A' . (5.9) 

Then, for j = 0, . . . , g, we have \\A^. |1 > and 

\Pnk, I ^ I dct Ak^ I - (C„,d,„fe, - l)v\\Ak, rSj . 

Proof. Let fc = 0, . . . , d. Using the multilinearity of the determinant we get 

p(A) ^det^^A.A^- =^det(AW),yl22)'---'^i(i))^^"^^'*'"^ ' 

where the sum is taken over all maps (p : {1, . . . , n} — > {0, . . . ,d}. Denoting by <i>fe 
the set of all such maps that satisfy J2m=i 4>{fn) — "^j we obtain that the nfc-th 
coefficient of the polynomial p is equal to: 

= detA,+ Y. det(A«),Ag)2),...,Aji)) . 

Using Hadamard's inequality together with the definition of rj yields 

I det(<(l), <(!,), . . . , a£))| ^ \\A%^h ■ ■ ■ Il<i)ll2 ^ ^P,(i)|| . . . P,(„)|l . 
Assume now that k = kj for some j = 0, . . . , q. Then, k is the abscissa of a vertex of 



the N ewto n polygon of the tropical polynomial Xp defined in Theorem 2.1 By Propo 
sition 



3.3 



applied to tp, we get that \\Au\\ > 0, \\A^\\ < ||Afc|| JlLm+i "<? ^ ll^fclkfe" 

for ah m ^ fc and \\A^\\ ^ \\Ak\\ Utk+i «7' < \\M4+T < Pfclkfc"'^?"" ^^r all 
m > k. Hence 

11^0(1)11 ■••||^0(«)ll ^ ll^fell ttfe Oj 

^ p^ll„^i:™=i,0(,„)>fc(</'M-fc) 

when X]m=i '/'('7^) — nk. When in addition <j) ^ k, there exists m — 1, . . . ,n such that 
(t){m) > k, thus Em=i,0(m)>fc('^("^) - fc) ^ 1: which yields 

p^(i)ii...iiA^(„)KPfeir5, . 

From all the above inequalities, we deduce that 

\Pnk\ ^ IdetAfel- Y, vWAkrSj ■ 

Since by definition Cn,d,nk is the cardinality of $fc, and there exists exactly one element 
of $fe such that (j) = k, we obtain the inequality of the lemma for k ^ kj. D 
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Proof of Theorem 2.12 Consider the polynomial p of (5.9) and let 71, . . . ,7nd 

maxo^;^„d \pi\x^ arranged in 



be the tropical roots of the tropical polynomial tp{x 
non-decreasing order. Let j = 0, . . . , 9 and denote k 



kj. From (1.2 1, we have 



ICl---CnfcK („",J71 



nd\ 
nk, 



■ • 7nfc 



(5.10) 



By the first part of Proposition 3.3 applied to the tropical polynomial tp, we get 
|p„fc|7i . . . 7nfc ^ |i5o|. Applying the result of Lemma 5.8 and using the assumption of 

\Ak\ 



Theorem |2. 12 
that 7i . . . 7„fc 



on Cj , we get that Pnk ^ Cj 



> 0. Since po = \ det Ao|, we deduce 



< 



|detAo|/(cjPfc||"'). Gathering this with (5.101, we get 



fnd\ 



ICl ■ • -Cnfel ^ 



IdetAolllAoll" 
Poll" P,||" 



Since k = kj is the abscissa of a vertex of the Newton polygon of the tropical polyno- 
mial tp, we get from Proposition 



3.3 



previous inequality shows Theorem 2 



applied to tp that 



Mo 

IIAfc 



= «! . . . a^, hence the 



12 D 



6. Examples. In this section we provide some examples to show that the lower 
bounds stated in Section [2] can be tight when the tropical roots are well separated and 
when the input matrices are well conditioned. We used the tropical scaling algorithm 
introduced in |GS09j to compute the eigenvalues of a given matrix polynomial. We 
warn the reader that a naive double precision numerical computation (linearization, 
followed by QZ algorithm) fails to give accurate values for some eigenvalues, because 
eigenvalues have different orders of magnitude (hence, either adapted scalings, or 
extended precision arithmetic, must be used). All the computations were performed 
in Scilab 5.3.0. 



In the following tables, we are comparing the value of the ratio 



ICi-.-C^fcl 



with 



(ai-.-ai-)" 

the lower bounds (2.1 1 of Section [2] for some matrix polynomials P. In view of Re- 
mark 5.7 we computed only the lower bounds obtained under (A2) for the normalized 
Frobenius norm || • ||„f, and the lower bounds obtained under |(Al)| for the normal- 
ized Schatten 1-norm (or normalized trace norm) || • \\^,i (which satisfies |(Al)| but not 
|(A2)[ see Property 5.4). Then, in each table, for each possible value of A; = 1, . . . , d, 
the column RATIO gives } \\i , and the other columns give the lower bound cVi 

of (2.1) with different values of Lk- The column G-p6lya gives the universal lower 



bound, generalizing the lower bound of Polya to the matrix case, for which Lk is given 



in (2.5) or (2.6), the column fewnom gives the lower bound involving the number 
of monomials, for which Lk is given in (2.2) or (2.3), and the column SEPAR gives 
the lower bound based on the separation between the tropical roots, for which Lj, is 



replaced by the constant L^'^°^ of (2.18) or 2^*'?''°'^ of (2.19) 



Example 6.1. Consider the following matrix polynomial 



PiiX)^10-'Uo 



lO^A^C/a 



lO^A'^C/4 + lOA^C/7 



lO-^A^C/g 



where all the matrices Uj, j E {0,2,4,7,9}, are unitary of dimension 3. Here, we 
shall only consider the normalized Frobenius norm since for a unitary matrix U, we 
have ||C/||nF = ||C^||*i = Ij so the tropical polynomial associated to Pi are the same 
for both norms, and the lower bounds under [(Al)| are the same, and are thus weaker 
than the lower bounds under |(A2)| for the normalized Frobenius norm. The Newton 
polygon of the tropical polynomial corresponding to Pi (for the normalized Frobenius 
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Fig. 6.1. The Newton polygon of the tropical polynomial corresponding to Pi. 



norm) is shown in Figure 6.1 and its tropical roots are 10 ^/^, 10 ^/^, 10^, and 10^^^ 
with multiphcities 2, 2, 3, and 2, respectively. Table [O] shows the lower bounds for 



all values of 1 ^ A: ^ 9, and the maximum of the ratios for a sample of 1000 random 
choices of the unitary matrices Uj. The ratios slightly change for different random 
choices of the unitary matrices Uj, but with a difference smaller than 5.910""'. Note 
that when 1 ^ k ^ 9 increases, the generalized lower bound of Polya decreases, and 
should be replaced by its symmetrized version with k replaced by d— k. However the 
bounds using the separation between tropical roots are the best ones, and they are 
tight at the vertices of the Newton polygon of the tropical polynomial associated to 
Pi. 



k 


RATIO 


FEWNOM 


G-POLYA 


SEPAR 


1 


1.00015 


0.08944 


0.125 


0.125 


2* 


1.00029 




0.05702 


0.97044 


3 


1.00015 




0.03425 


0.125 


4* 


1 




0.02345 


0.99991 


5 


1.00000 




0.01733 


0.125 


6 


1.00000 




0.01348 


0.12500 


7* 


1 




0.01088 


0.99056 


8 


1 




0.00901 


0.125 


9* 


1 




0.00763 


1 



Table 6.1 
Comparison of the bounds for Pi using the normalized Frobenius norm. The abscissiE of the 
vertices of the Newton polygon of corresponding tropical polynomial are indicated by the symbol * . 



Example 6.2. In the following example, we perturb the previous polynomial 
by adding coefficients to Pi , in such a way that the Newton polygon of the tropical 
polynomial associated to Pi does not change, either for the normalized Frobenius norm 
or at least for the normalized Schatten 1-norm. Here we fix the unitary matrices to 
be either the identity matrix or its opposite, and wc consider the following matrices: 



B^b 



10-1^ 


1 


1 





1 


1 


1 


2 


3 



c = 




where b is chosen so that ||P||nF = l, and C is also such that ||C||nF=l- The matrix 
C is singular, since it has rank 1, and it satisfies ||C||*i — l/\/3- The matrix B is 
ill-conditioned, so that llPlLi < ||P||nF == 1 too. 
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We shall consider the two polynomials: 

P2(A) ^IQ-^Uo + IQ-^/^ACi + 1Q^X^U2 + lO^/^X^Ci + lO'^X'^U^ + lO^A^Ci 

+ lO^A'^Bi + lOA^t/7 + IQ-'^/^X^Bi + IQ-^X^Ug 
P3(A) ^IQ-^Uo + IQ-^/^ACa + lO^A^t/a + 10^/^X^02 + lO'^X'^Ui + lO^A^Cz 

+ 10^X^B2 + lOA^t/y + lO^^/^A^Sa + IQ-^X^Ug , 



with Uo ^ U4 = U7 = -I, U2 = Uq = I, Bi ^ B, Ci = C, B2 = \\B\\ 



'B, 



C2 = vSC*. Since for any unitary matrix U with dimension 3, ||f7||*i = ||f^||nF — Ij 
and since for any complex matrix A of dimension 3, we have ||^||*i ^ ||A||nF, the 
tropical polynomial associated to P2 for either the normalized Frobenius norm or 
the Schatten 1-norm coincides with the one associated to Pi. Hence, as for Pi, the 
bounds for P2 with the normalized Schatten 1-norm are necessarily weaker than the 
ones with the normalized Frobenius norm. Moreover, all the bounds presented in 



Table 6.1 for Pi remain the same for P2 except the bounds based on the number of 



nonzero coefficients. We present all these bounds in Table 6.2 together with the new 
ratios. Since these ratios are different from the ones of Table |6.1[ one can see that 
the lower bounds based on the separation between tropical roots may be tight. 



k 


RATIO 


FEWNOM 


G-POLYA 


SEPAR 


1 


0.45628 


0.03162 


0.125 


0.125 


2* 


1.03298 




0.05702 


0.97044 


3 


0.45627 




0.03425 


0.125 


4* 


1.00009 




0.02345 


0.99991 


5 


0.39258 




0.01733 


0.125 


6 


0.39258 




0.01348 


0.12500 


7* 


1.00867 




0.01088 


0.99056 


8 


0.47428 




0.00901 


0.125 


9* 


1 




0.00763 


1 



Table 6.2 

Comparison of the bounds for P2 using the normalized Frobenius norm. The abscissds of the 
vertices of the Newton polygon of the corresponding tropical polynomial are indicated by the symbol 



Let us consider now the polynomial P3. The normalized Schatten 1-norm of the 
coefficients of P3 coincide with the normalized Frobenius norm of the coefficients of 
P2 , hence the tropical polynomial associated to P3 for the normalized Schatten 1-norm 
coincides with the one associated to Pi for the normalized Frobenius norm, so that 
the tropical roots are the same. However, since the Schatten 1-norm satisfies |(Al)| 
but not |(A2)[ we can only get the bounds based on |(Al)| which are necessarily lower 
than the one presented above for Pi with the normalized Frobenius norm. Finally, the 
tropical polynomial associated to P3 for the normalized Frobenius norm differs from 
the one associated to Pi, so that the tropical roots and the bounds differ too. All 
tropical roots of this new tropical polynomial have multiplicity 1, which means that 
all indices are abscissa of vertices of its Newton polygon. We present in Table |6.3[ 
the ratios and the lower bounds for both normalized Frobenius and Schatten 1-norms. 
Here the ratios and lower bounds for the normalized Schatten 1-norm are lower than 
the ones for the normalized Frobenius norm. However, one can see that the results 
for the normalized Frobenius norm are still better than the ones for the normalized 
Schatten 1-norm. 

Example 6.3. In order to obtain a polynomial for which the lower bounds for 
the normalized Schatten 1-norm are better than the ones for the normalized Frobenius 



26 



M. Akian , S. Gaubert and M. Sharify 





norm. Frobenius norm 


norm. Schatten 1-norm 


k 


RATIO 


FEWNOM 


G-POLYA 


SEPAR 


RATIO 


FEWNOM 


G-POLYA 


SEPAR 


1 


1.57223 


0.03162 


0.125 


0.35355 


0.30258 


0.001 


0.01563 


0.01563 


2* 


1.11466 




0.05702 


0.91391 


1.11466 




0.00325 


0.54771 


3 


1.57223 


" 


0.03425 


0.35355 


0.30258 


" 


0.00117 


0.01562 


4* 


1.00028 




0.02345 


0.99972 


1.00028 




0.00055 


0.96682 


5 


1.45041 


" 


0.01733 


0.18345 


0.27913 


" 


0.00030 


0.01562 


6 


1.01082 


" 


0.01348 


0.05989 


0.27913 


" 


0.00018 


0.01533 


7 


1.02104 




0.01088 


0.97788 


1.02104 




0.00012 


0.71337 


8 


1.26686 


" 


0.00901 


0.25721 


0.34983 


" 


0.00008 


0.01563 


9* 


1 




0.00763 


1 


1 




0.00006 


1 



Table 6.3 
Comparison of the bounds for P3 using the normalized Frobenius and 
abscissas of the vertices of the Newton polygon of the tropical polynomial 
normalized Frobenius norm are indicated by the symbol *. 



Schatten 1-norms. The 
associated to P3 for the 



norm, one need to take matrices Ai such that ||Ai|jnF/||^i||*i is as large as possible. 
This means that the dimension n is large and that some of the singular values of A^ 
are zero, or at least that the ratio between the maximal and minimal singular values 
is large, which implies that Ai is singular or at least ill-conditioned. Moreover some 
parameters of the lower bounds need to be smaller than ||Ai||nF/||^i||*i, which is 
possible for instance for the lower bounds based on the separation between tropical 
roots, when these tropical roots have a small multiplicity and arc well separated. We 
shall show here such examples. 

We consider a matrix polynomial with degree 3 and dimension n = 10: 

PiiX) = Ao + AiA + A2X^ + AsX^ , 

where all the matrices Ai,i ~ 0, . . . , 3 are diagonal matrices: Ai ~ diag((ii) with 
di e C". We choose Aq = A3 = /, so that c = 1, a\(X20i-i — 1 (for norms such that 
II /|| = 1) and for A: = 3 the ratio }^" ,^^\n is equal to 1, which is equal to the lower 
bound based on separation of tropical roots. We choose the diagonal elements of the 
matrices A 1,^2 as follows: 



with ai,a2,b > 0. Then 



di:= 
d2:= 



ai 



1 + b/n 

Q2 

1 + b/n 



(l,...,l,l + fe) , 
(1 + 5,1,...,1) , 



4 11 ll^llll 
|Ai||*i = = ai 



\A. 



^2||l 



= 02 



and 



1 A II ll'^llb r 

\Ai\\nF = 1=^ = aijb 



\A 



2||nF 



12\\2 



= 0-2 fb 



with 



Jb ■= 



2b+b^ 



l + b/n 



> 1 
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The constant fb is large when b is large, that is when the condition numbers of Ai 
and A2 are large. 

When oi = a2 > 1, the tropical polynomial associated to P4 for the normalized 
Schatten 1-norm (resp. the normalized Frobenius norm) has 3 different tropical roots 
equal to l/oi, 1 and oi (resp. l/(ai/f,), 1 and ai/b). We present in Tables 6.4 and 6.5 



the ratios and the lower bounds for both normalized Frobenius and Schatten 1-norms, 
when fli — 10'', and 6 = 2 and 6 = 10 respectively. When b — 10, the matrices are 
ill-conditioned, so that the ratios J "';L for fc = 1 or 2 are very large. When either 
& = 2 or & = 10, we see that the lower bounds based on the separation between the 
tropical roots are the best ones, and that they are nearer the ratios ( J^^'"^"^„ ) in the 
case of the Schatten 1-norm than in the case of the Frobenius norm, although they 
are far in the case where 6 = 10. 





norm. Frobenius norm 


norm. Schatten 1-norm | 


k 


RATIO 


FEWNOM 


G-POLYA 


SEPAR 


RATIO 


FEWNOM 


G-POLYA 


SEPAR 


1* 
2* 
3* 


6.30697 
6.30697 

1 


0.00098 
0.00098 
0.00098 


0.00098 
0.00007 
0.00001 


0.99911 
0.99911 
1 


2.06667 
2.06667 
1 


9.5D-07 
9.5D-07 
9.5D-07 


9.5D-07 
5.1D-09 
1.7D-10 


0.81873 
0.81873 
1 



Table 6.4 
Comparison of the bounds for P4 with ai = 02 = 10* and 6 = 2. The abscissas of the vertices 
of the Newton polygon of corresponding tropical polynomial are indicated by the symbol * . 



norm. Frobenius norm 

RATIO H FEWNOM \ G-POLYA 



RATIO 



norm. Schatten 1-norm 

FEWNOM \ G-POLYA I SEPAR 



1* II 33882.7 II 0.00098 I 0.00098 I 0.99945 11 93.4462 11 9.5D-07 I 9.5D-07 0.81873 
2* 33882.7 0.00098 0.00007 0.99945 93.4462 9.5D-07 5.1D-09 0.81873 

3* II 1 II 0.00098 I 0.00001 | 1 |_JL | 9.5D-07 | 1.7D-10 | 1 

Table 6.5 
Comparison of the bounds for P4 with ai = 02 = 10'' and b = 10. The abscissas of the vertices 
of the Newton polygon of corresponding tropical polynomial are indicated by the symbol * . 



When tti > 1 and 02 = \fa\ffh the tropical polynomial associated to P4 for 
the normalized Schatten 1-norm (resp. the normalized Frobenius norm) has only 2 
different tropical roots equal to l/oi and ^Jli{ (resp. l/(ai/b), 1 and \faiYh) with 



respective multiplicities 1 and 2. We present in Tables 6.6 and |6.7[ the ratios and the 
lower bounds for both normalized Frobenius and Schatten 1-norms, when a\ — 10^, 
and 6 = 2 and 6 = 10 respectively. When 6 = 10, the matrices are ill-conditioned, so 
that for fc = 1 or 2 the distances (ratios) between the ratio }}-^ and its lower 
bounds are very large, although for A: = 2 this ratio is smaller than in the case where 
a\ — a2 — 10^ above. When either 6 = 2 or 6 = 10, we see that the lower bounds 
based on the separation between the tropical roots are the best ones. For k — 2 which 
is the abscissa of a vertex of the Newton polygon, the lower bound obtained in the 
case of the Schatten 1-norm is nearer the ratio , ^„ than the one obtained in the 

(qi. ..«/,)" 

case of the Frobenius norm, whereas for fc = 3 the contrary holds. 
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norm. Frobenius norm 
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